Entanglement as Measure of Electron-Electron Correlation in 
Quantum Chemistry Calculations 

Zhen Huang and Sabre Kais*^ 
Department of Chemistry, Purdue University, West Lafayette, IN 4^907 

Abstract 

In quantum chemistry calculations, the correlation energy is defined as the difference between 
the Hartree-Fock limit energy and the exact solution of the nonrelativistic Schrodinger equation. 
With this definition, the electron correlation effects are not directly observable. In this report, 
we show that the entanglement can be used as an alternative measure of the electron correlation 
in quantum chemistry calculations. Entanglement is directly observable and it is one of the most 
striking properties of quantum mechanics. As an example we calculate the entanglement for He 
atom and H2 molecule with different basis sets. 
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The Hartree-Fock self-consistent field approximation, which is based on the idea that we 
can approximately describe an interacting fermion system in terms of an effective single- 
particle model, remains the starting point and the major approach for quantitative elec- 
tronic structure calculations. In quantum chemistry calculations, the correlation energy is 
defined as the energy error of the Hartree-Fock wave function, i.e., the difference between the 
Hartree-Fock limit energy and the exact solution of the nonrelativistic Schrodinger equationi. 
There also exists other measures of electron correlation in the literature such as the statistical 
correlation coefficients^ and more recently the Shannon entropy as a measure of the correla- 
tion strength-'^. Electron correlations have a strong influence on many atomic, moleculai?^, 
and solid properties^. Recovering the correlation energy for large systems remains one of 
the most challenging problems in quantum chemistry. 

The concept of electron correlation as defined in quantum chemistry calculations is useful 
but not directly observable, i.e., there is no operator in quantum mechanics that its mea- 
surement gives the correlation energy. In this letter, we propose to use the entanglement as 
a measure of the electron correlation. Entanglement is directly observable and it is one of 
the most striking properties of quantum mechanics. 

It was nearly 70 years ago when Schrodinger gave the name "entanglement" to a correla- 
tion of quantum nature. He stated that for an entangled state "the best possible knowledge 
of the whole does not include the best possible knowledge of its parts"''. Over the decades 
the meaning of the word "entanglement" has changed its flavor and our view of the na- 
ture of entanglement may continue to be modified^. Entanglement is a quantum mechanical 
property that describes a correlation between quantum mechanical systems that has no clas- 
sical analog^iiSiiiii^. A pure state of a pair of quantum systems is called entangled if it is 
unfactoriazable, as for example, the singlet state of two spin-^ particles, -^(1 t|) — | it))- 
A mixed state is entangled if it can not be represented as a mixture of factorizable pure 
statesi^ii^^. Since the seminal work of Einstein, Podolsky, and Rosei>i^ there has been a 
quest for generating entanglement between quantum particlesiSiii. Investigation of quan- 
tum entanglement is currently a very active area and has been studied intensely due to 
its potential applications in quantum communications and information processings^ such as 
quantum teleportatio n^^i^^ , superdense coding^°, quantum key distribution^^, telecoloning, 
and decohernce in quantum computers^^'^S'^i. 

In order to employ the entanglement as an alternative method of measuring the electron 
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correlation, we present the scheme to quantify entanglement based on the entanglement 
measure for two-particle system o^^i^^ . We will obtain a general approach to quantify the 
entanglement between different spin-orbitals of atomic and molecular systems. 

For two electron system in 2m-dimensional spin-space orbital with and denote the 
fermionic annihilation and creation operators of single-particle states and |0 > represents 
the vacuum state, a pure two-electron state |$ > can be written as 



(1) 



a,6e{l, 2,3,4,. ..2m} 

where a, b run over the orthonormlized single particle states, and Pauli exclusion requires 



that the 2mx 2m expansion coefficient matrix u is antisymmetric: u, 



a,b 



-ujbn, andcuji = . 



In the occupation number representation (ni 1,^2 I,..., rim h^m I), 

where t and | mean a and (3 electrons respectively, the subscripts denote the spatial orbital 
index and m is the total spatial orbital number. By tracing out all other spatial orbitals 
except rii, we can get a (4x4) reduced density matrix for the spatial orbital rii 
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where Pmfi denotes an "empty orbital" , 

m—l m—1 

Pnifi = 4 E E l^2i+l,2i+2| 
i=l j=l 

and pni,2 denotes "two electron occupied orbital", 

Pnl,2 = M^lM^- 



(2) 



(3) 



(4) 



(5) 



The "one electron occupied orbital" , in (|, |) basis set is given by 

4EIllM^2,2ml' 

^ET=2M 

The matrix elements of u can be calculated from the expansion coefficient of the ab initio 
Configure Interaction (CI) method. The CI wave function with single and double excitation 
can be written as 
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ar a<b,r<s 

where |\l/o > is the ground state Hartree-Fock wave function, cj^ is the coefficient for single 
excitation from orbital a to r, and c^'^ is the double excitation from orbital a and 6 to r and 
s. Now the matrix elements of uj can be written in terms of the CI expansion coefficients 

J2i+1 2i+2 2i+l,2j+2 
_ Co _ Ci _ C2 _ C12 , . 

^1,2 — —, t^2,2i+l — 2 ' ^l,2j+2 — 2 ' '^2i+l,2i+2 — ' , [1 ) 

where i,j = l,2...m. In this general approach, the ground state entanglement is given 
by von Neumann entropy of the reduced density matrix p„i 

S{prn) = -Tr{prJog2prn)- (8) 

We are now ready to evaluate the entanglement for the H2 molecule as a function of R 
using a direct and simpler approach based on the two-electron density matrix calculated 
from the CI wave function with single and double electronic excitations. In the occupation 
number representation, the CI wave function is given by 



1$ >= colllOO... > +c?|0110... > +c^|1001... > +c?;2|0011... > +.... (9) 

By tracing out all other orbitals except 1, we can get the reduced density matrix for {rii |= 
0,1) 

(Y^m-l |^2j+l|2 I i 2i+l,2i+2|2 r, \ 

The CI wave function expansion coefficients are calculated with the electronic structure 
package GaussianSi. Thus the entanglement of H2 molecule is readily calculated by the von 
Neumann entropy 

Sipi) = -Tr{pilog2pi). (11) 

Figure (1) shows the calculated entanglement S for II2 molecule, Eq. (jllll . as a function 
of the internuclear distance R using Gaussian basis set 3-21 G^^. For comparison we included 
the usual electron correlation {Ec = \E^^°''^^ — E^^^\) and spin- unrestricted Hartree-Fock 
(UHF) calculations^! using the same basis set in the figure. At the limit R = 0, the dot 
represents the electron correlation for the He atom, Ec = 0.0149(a.u.) using 3-21G basis 



set compared with the entanglement for the He atom 5* = 0.0313. With a larger basis set, 
cc — pV5Z'^^, we obtain numerically Ec = 0.0415(a.M.) and S = 0.0675. Thus, qualitatively 
entanglement and absolute correlation have similar behavior. At the united atom limit, 
R 0, both have small values, then rise to a maximum value and finally vanishes at the 
separated atom limit, i? — > oo. However, note that for R > 3 A the correlation between the 
two electrons is almost zero but the entanglement is maximal until around R ~ 4A, the 
entanglement vanishes for R > 4A. 

To understand the entanglement behavior for H2 molecule using ab initio quantum chem- 
istry methods, we calculate the entanglement for a simpler two-electron model system. This 
is a model of two spin-i electrons with an exchange coupling constant J (a.u.) in an effective 
transverse magnetic field of strength B (a.u.). In order to describe the environment of the 
electrons in a molecule, we simply introduce a small effective external magnetic field B. The 
general Hamiltonian for such a system is given by 



H = -^{1 + -f)a^ 0a^- ^{1 - 7)af ® - Ba^ h - Bh ® (^2, 



(12) 



where a" are the Pauli matrices(a= x,y,z) and 7 is the degree of anisotropy. For 7 = 1 
Eq. lfT^ reduces to the Ising model, whereas for 7 = it is the XY model. 

Our two spin problem admits an exact solution, it is simply a (4 x 4) matrix with the 
following four eigenvalues 



Ai = -J, A2 = J, A3 = -^452 + J272, A4 = + J272 

and the corresponding eigenvectors 
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(14) 



where a = a/4B^~+"J^7^. In the basis set {| tT>) I Ti>) I iT>) I ii>}) the eigenvectors can 
be written as 



H >=_(iit>+in>) 



(15) 



/>2>=4^(iiT>-in>), (16) 



a-2B. , , a + 2B. 



Now we confine our interest to tlie calculation of tlie entanglement between tlie two 
electronic spins. For simplicity we take 7 = 1, Eq. ()12|) reduces to the Ising model with 
the ground state energy A3 and the corresponding eigenvector |03 >. All the information 
needed for quantifying the entanglement in this case is contained in the two-electron density 
matrix. 

When a biparticle quantum system AB is in a pure state there is essentially a unique 
measure of the entanglement between the subsystems A and B given by the von Neumann 
entropy S"^^. If we denote pa the partial trace of pab with respect to subsystem B, pA = 
TrsipAB), the entanglement of the state pab is defined as the von Neumann entropy of the 
reduced density operator pA, S{pab) = —Tr[pAlog2PA]- 

For our model system in the ground state \4>3 >, the reduced density matrix in the basis 
set(t, i) is given by 



/ a+2h 



Pa 



2a 

0a—2h 



(19) 



Thus, the entanglement is simply given by 



where X = J/B. 

The value of J, the exchange coupling constant between the spins of the two electrons, 
can be calculated as half the energy difference between the lowest singlet and triplet states 
of the hydrogen molecule. Herring and Flicker have shown^ that J for H2 molecule can be 
approximated as a function of the interatomic distance R. In atomic units, the expression 
for large R is given by 

J{R) = -0.821 i?^/2g-2K ^ o(^2g-2fl^)_ (21) 

Figure (2) shows the calculated von Neumann Entanglement (S), Eq. ()20|) . as a function 
of the distance between the two electronic spins R, using J{R) of Eq. (pT|) . for different 



values of the magnetic field strength B. At the limit i? — oo the exchange interaction J 
vanishes as a result the two electronic spins are up and the wave function is factorazable, 
i.e. the entanglement is zero. At the other limit, when R = the entanglement is zero 
for this model because J = 0. As i? increases, the exchange interaction increases leading 
to increasing entanglement between the two electronic spins. However this increase in 
the entanglement reaches a maximum limit as shown in the figure. For large distance, 
the exchange interaction decreases exponentially with R and thus the decrease of the 
entanglement. The figure also shows that the entanglement increases with decreasing the 
magnetic field strength. This can be attributed to effectively increasing the exchange 
interaction. Thus we get the similar behavior as the entanglement for the H2 molecule as a 
function of the internuclear distance R , using accurate ab initio methods. Because the Eq. 
(PT|) is only applicable to the large value of R, it is not surprising to see, in the hmit — > 0, 
the entanglement converges to the He atom results in ab initio methods but disappears in 
this simple model. 

Recently a new promising approach is emerging for the realization of quantum chemistry 
calculations without wave functions through first order semidefinite programming^. The 
electronic energies and properties of atoms and molecules are computable simply from an 
effective two-electron reduced density matrix P(ab)- Thus, the electron correlation can be 
directly calculated as effectively with the entanglement between the two electrons, which is 
readily calculated as the von Neumann entropy S = —TrpAlog2PA, where pA = TrBP(AB)- 
Utilizing this combined approach, one calculates the electronic energies and properties of 
atoms and molecules including correlation without wave functions or Hartree-Fock reference 
systems. 

In summary, we presented the entanglement as an alternative measure of the electron 
electron correlation in quantum chemistry calculations for atoms and molecules. All the 
information needed for quantifying the entanglement is contained in the two-electron density 
matrix. This measure is readily calculated by evaluating the von Neumann entropy of the one 
electron reduced density operator. This definition of correlation has deep roots in quantum 
theory, observable, and does need a reference system such as Hartree-Fock calculations. The 
approach is general and can be used for larger atomic and molecular systems. 
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FIG. 1: Comparison between the absolute value of the electron correlation Ec = \E^^^'^* — 
E^^^\ and the von Neumann Entanglement (S) as a function of the internuclear distance R for 
the H2 molecule using Gaussian basis set 3-21G. At the limit R = 0, the dot represents the 
electron correlation for the He atom, Ec = 0.0149(a.it.) using 3-21G basis set compared with 
the entanglement for He atom S = 0.0313. The equilibrium distance using 3-21G basis set is 

Ren = O.74A 
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FIG. 2: von Neumann Entanglement (S) as a function of the distance R between the two spins for 
different values of the magnetic field strength B. 
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